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Abstract 

The impact of environmental variability on population size growth rate in dynamic 
models is a recurrent issue in the theoretical ecology literature. In the scalar case, 
R. Lande pointed out that results are ambiguous depending on whether the noise is 
added at arithmetic or logarithmic scale, while the matrix case has been investigated 
by S. Tuljapurkar. Our contribution consists first in introducing another notion of vari- 
ability than the widely used variance or coefficient of variation, namely the so-called 
convex orders. Second, in population dynamics matrix models, we focus on how matrix 
components depend functionaly on uncertain environmental factors. In the log-convex 
case, we show that, in a sense, environmental variability increases both mean popu- 
lation size and mean log-population size and makes them more variable. Our main 
result is that specific analytical dependence coupled with appropriate notion of vari- 
ability lead to wide generic results, valid for all times and not only asymptotically, and 
requiring no assumptions of stationarity, of normality, of independency, etc. Though 
the approach is different, our conclusions are consistent with previous results in the lit- 
erature. However, they make it clear that the analytical dependence on environmental 
factors cannot be overlooked when trying to tackle the influence of variability. 

Key words: environmental variability; matrix population models; growth rate; 
stochastic orders; log-convex functions. 
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1 Influence of environmental noise on population size 

We recall here different observations and results in the theoretical ecology literature which 
point out the ambiguous role of environmental noise on population size in matrix population 
models, according to whether the noise is added at arithmetic or logarithmic scale. 

1.1 Lande's comments on additive noise at arithmetic or logarith- 
mic scale 

R. Lande in [Ij comments the influence of environmental noise on population size according 
to whether the noise is added at arithmetic or logarithmic scale. The evolution of population 
size N{t) in absence of density-dependent effect may be described 

• either on arithmetic scale with multiplicative growth rate X{t) and dynamic 
N{t + 1) = Xit)N{t), 

• or on logarithmic scale with growth rate on the log scale r(t) = log (A(t)) and dynamic 
on the log scale log A^(t + 1) = r{t) + log N(t). 

On the one hand, adding environmental noise to multiplicative growth rate as in A(t) = 
A + e{t), where the noise is zero-mean (E[e(t)] = 0), gives the following mean of growth rate 
on the log scale 

r = E[log (A(t))] ^logA-a^ . 

"Thus, demographic and environmental stochasticity reduce the mean growth rate of a pop- 
ulation on the logarithmic scale, compared with that in the (constant) average environment" 
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On the other hand, adding environmental noise to growth rate on the log scale as in 
r{t) = f + e(t) gives, in case e(t) follows a Normal distribution A/'(e, a^), the following mean 
of growth rate on the arithmetic scale 

A = exp (f + e + -;^) . 

Thus, Lande concludes that, "with the mean environmental effect equal to zero, e = 0, then 
it would be found that environmental stochasticity increases the mean multiplicative growth 
rate. A". 



1.2 Tuljapurkar's asympotic approximation 

S. Tuljapurkar considers a stationary sequence of random matrices Aq, Ai, . . .yielding pop- 
ulation vector n{t) = At-\ ■ ■ ■ Aonifi) and population size N{t) = \\At-i ■ ■ ■ y4on(0)||. Under 
general conditions (see [21 [3]), there exists a deterministic stochastic growth rate As defined 
by 

logA,= lim -logN{t) = lim - log || A_i • ■ ■ Aon(0)|| . 

t^+oo t t— >+c« t 

Denoting by Ai the largest eigenvalue of the average matrix A, Tuljapurkar obtains the 
approximation 



logA, ^logAi-— 2 + ^ 



9_ 

2Af ^ Af 

where is proportional to the variance E[(y4t — A) (At — A)] (and 9 is related to auto- 
correlation). In this case, environmental stochasticity reduces the mean growth rate of the 
population. 



1.3 A quest for generic results 

The two above cases show that environmental noise has an ambiguous impact on population 
size in matrix population models. Our main objective is contributing to clarify this impact 
with generic mathematical results. For this, we shall first introduce in Sect. [2] a tool to mea- 
sure variability, distinct from the widely used variance or coefficient of variation, and known 
as convex partial orders. Then, in Sect. El we shall provide generic results on environmental 
noise variability in population dynamics matrix models. We conclude in Sect. Hlby pointing 
out proximities and differences between our approach and those presented in Sect. |2l 



2 Convex orders as tools for measuring variability 

To a (square integrable) random variable X, one can attach the variance var(X). This latter 
scalar measures "variability" , and any pair of random variables X and Y may be compared, 
with X being more variable than Y if var(X) > var(y). The variance thus defines a total 
order. 
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Other orders are interesting for comparing pairs of random variables. However, they are 
generally not total: not all pairs may be ranked. Related to this is the fact that no single 
scalar, such as variance, may be attached to a random variable to measure its variability. In 
this vein, we shall present the so-called increasing convex and convex stochastic orders. Such 
orders can only rank random variables for which the primitives of their respective repartition 
functions never cross. 

We think that these orders and many others referenced in the two main books [U [S] may 
be useful in the ecological modelling scientific community. Of course, for this, practical tests 
must be developed to compare empirical data as to their variability. This is not the object 
of this paper. 

All random variables are defined on a probability space with probability P. To a random 
variable X, we shall attach its (right-continuous) repartition function F{x) = P(X < x). 
We shall always consider random variables with finite means, with generic notation X and 
Y, and F and G for their respective repartition functions. 

2.1 Increasing convex order 

The increasing convex order compares random variables according both to their "location" 
and to their "variability" or "spread" [5]. We say that X is less than Y in increasing convex 
order, denoted by 

Y ^ V 

^ —icx -' ) 

if and only if one of the following equivalent conditions holds true 

• the primitive of the repartition function of X is always below that of Y: 
/f ^ F{x)dx < ^ G{x)dx, for all c G M, 

• E((/?(X)) < E((p(y)) for all increasing and convex function cp. 

Roughly speaking, Y is more likely to take on extreme values than X. In a sense, X is 
both "smaller" and "less variable" than Y |5]. We have the important property that, when 
^ :^icx Y 1 the means are ordered too: E(X) < E(y). However, nothing can be said of the 
variances. To compare variances, we need a stronger (more demanding) order. 

2.2 Convex order 

The convex order compares random variables according to their "variability" or "spread" 
[5]. We say that X is less than Y in convex order, denoted 

X :<cx Y , 

if and only if one of the following equivalent conditions holds true: 

• the means are equal and the primitive of the repartition function of X is always below 
that of Y, that is, E(X) = E(F) and X Y, 
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(a) Repartition functions 



(b) Primitives 



Figure 1: Illustration of the convex order: the primitives of the repartition functions do not 
cross each other 

• E,{ip{X)) < E,{ip{Y)) for all convex function (p. 

Roughly speaking, Y is more likely to take on extreme values than X (see Figured]). Notice 
that the convex order is more demanding than the increasing convex order since the class 
of "test functions" is larger: all convex functions and not only the increasing convex ones. 
This is why we obtain stronger important properties that, when X :<cx Y, the means are 
equal E(X) = ]E(y), and the variance are ordered var(X) < var(y). 

2.3 Some properties 

• This icx and cx orders are stricter than the order defined by comparing variances: not 
all pairs of random variables may be ranked. 

• Consider the class of random variables having same mean /x and variance cr^. 
Elements of 0-2 cannot be compared with respect to cx. Indeed, if X :<cx Y and 
var(X) = var(y), then X and Y have the same distribution [T, p. 57]. 

• Adding zero mean independent noise to a random variable increases variability: if Z 
is independent of X and has zero mean, then X is less than Y = X + Z in convex 
order. This is a consequence of Strassen's Theorem [ij p. 23]. More generally, without 
assuming independence, X is less than Y = X + Z in convex order whenever the 
conditional expectation E[Z|X] = 0. 

• Consider X following Normal distribution A/'(/i, o"^) and Y following M{v,t'^). Then, 
X :<icx Y if and only if /x < z/ and o"^ < r^, and X Y if and only if /i = z/ and 

< [4, p.62]. 

• For p > 0, let us introduce CVp(X) := ^(^x) positive p-integrable random vari- 
able X. For p = 2, we have the usual coefficient of variation CV(X) = CV2(X) = 
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v/E(X2)/E(X). If X then CV2(X) < CV2(F) (in fact CVp(X) < CVp(r) for 

all p > 1). 

2.4 Increasing convex order and convex order for random vectors 

We shall need to compare not only random variables but random vectors as in [H p. 98] and 
[5l p. 323]. For this, we can no longer appeal to repartition functions. Let X = (Xi, . . . , X„) 
and Y = (Yi, . . . , Yn) be random vectors with finite mean. 

We say that X is less than Y in increasing convex order, written X ■<icx Y , if and only 
if E(99(Xi, . . . , Xn)) < E(y9(Y'i, . . . , Yn)) for any increasing convex function if : R" — > R. 

We say that X is less than Y in convex order, written X -<cx Y , if and only if E(y9(Xi, . . . , X, 
E(93(Yi, . . . , for any convex function : R" R. In this case, X and Y have the same 
mean. 

Consider X following Normal distribution 7V(/i, S) and X' following A^(yu', E'). Then, 
X -<cx X' if and only if fi = fi' and S' — S is non-negative definite. The situation is not 
as clear cut for the icx order. If fix > A^y and Ex — > (non-negative definite), then 
X ^icx Y. If X ^icx Y, then fix > yUy and a"^(Sx — Sy)a > for all vector a > [4t p. 100]. 

3 Generic results on environmental noise variability in 
population dynamics matrix models 

In what follows, we shall consider a population described at discrete times t = 0, . . . ,T (where 
T is the horizon), either by a scalar n(t) G R or by a vector n{t) = (rii(t), . . . , nn{t)) E R" 
which may be abundances at ages or stages. The population size is N = \\n\\ = ni-\ — ■ + nn. 
The dynamical evolution of the population is supposed to be linear in the sense that 



where the matrix A is independent of n{t) (no density-dependence effect, this is why we label 
such model of linear). On the other hand, the components Aij of the matrix A may depend 



For instance, the components of the matrix A may depend linearly on the environmental 
factors, as in the expression 



n{t + 1) = A{e{t))n{t) , t = 0,...,T -1 



(1) 





or may depend exponentially as in 




(3) 
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In this latter case, the components of the matrix A are log-convex function of the environ- 
mental factors. Recall that / is a log-convex function if / > and log / is convex. Otherwise 
stated, / is the exponential of a convex function (as a consequence, a log-convex function is 
also convex). 

In |6], different non linear models are recalled. When a = 0, they are matrix models 
without density- dependency. Model (2a) exhibits components which are exponential in the 
environmental factor, while they are linear in models (2c) and (2d). Calculation shows that 
model (2b) has matrix components which are log-convex functions of the environmental 
factor. 

We shall coin environmental scenario a temporal sequence e{-) = (^(0), . . . ,e{T — 1)) of 
environmental factors. 

Proposition 1 Consider two environmental scenarii, one being more variable in increasing 
convex order than the other: (^^^(O), . . . , e^^{T — 1)) ^jcz (^^^(0), . • • , e^{T — 1)) . Denote by 
N^'{T) = \\A{e'\T-l))---A{e''\0))n{0)\\ and N'^iT) = || ^(^^(T - 1)) ■ ■ ■ A(£^(0))n(0)|| 
the corresponding populations sizes. 

Assume that the components Aij{ei, . . . ,ep) of the matrix A in are nonnegative com- 
binations of log-convex functions of the environmental factor (ei, . . . ^Sp). Then, the more 
variable the scenario, the more variable the population size in the sense that 

N^\T) N\T) and \ogN^\T) logiV^(T) . (4) 

As a consequence, E (A^*^(r)) > E {N^{T)) and E (logA^*^(T)) > E (logAr^(T)). 

In a sense, environmental variability increases both mean population size and mean log- 
population size and makes them more variable. 

Proof. The components of the vector n(r) =A{e{T-l)) • • • A(e(0))n(0) are sums of products 
of nonnegative combinations of log-convex functions of the environmental scenario. Therefore, by a 
property of log-convex functions [7], the components of the vector n(T) are also log-convex functions 
of the environmental scenario, and so is the population size. Thus, the logarithm log N(T) of the 
population size is convex in (e(0), . . . ,e(T — 1)). For any increasing convex function : M ^ M, 
ip[log N{T)) is convex in (e(0), . . . , e(T — 1)) since convexity is preserved by left-composition 
with an increasing convex function. We end up by using the definition of increasing convex order 
for random vectors in g231 E[ip{log N^^ {T))] > E[99(logiV'^(r))] . This precisely means that 
logiV^(r) hie. log N^{T). 

Since a log-convex function is also convex, the population size n(T) is a sum of convex functions 
of the environmental scenario. Then, the proof follows as above. 

At last, we use the property that X 'hicx Y =^ E(X) > E(y) to compare the means. 

□ 

Instead of total population, the result would still hold true with any positive weighted 
combination aiUi -|- ■ ■ ■ -|- akUk where a, > 0, or with log(aini -|- ■ ■ ■ -|- akUk) where > 0, 

As an illustration, consider the following scalar dynamic equation for population size 
n(t + l) = exp (^r + e(t))n(t) , for which we have n(T) = exp (^rT + e{0) + ■ ■ ■ + e(T — l)^n{0) . 
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Hence, both n(T) and \ogn(T) are convex functions of the environmental scenario e(-) = 
{e{0), . . . ,e{T — 1)), so that environmental variability increases mean population size as 
may be seen in Figure [H Indeed, the mean population size generated by a more variable 
environment is above the one by a less variable environment, for all times. 



Population in Logarithm Scaie (Exponentiai Modei) 
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Figure 2: Environmental variability increases mean population size for a population model 
where growth rate depends exponentialy on environmental factor 



4 Conclusion 

We have used another notion of variability than the widely used variance or coefficient of 
variation, namely the so-called convex orders. We think that such partial orders may be of 
interest in theoretical ecology beyond this specific application. 

To compare our approach with the literature, notice that, though we consider matrix pop- 
ulation models, we make no ergodic assumption on the stochastic process Aq, Ai, ... However, 
we make separate assumptions, on the one hand on the environmental factors ^(0), . . . , e{t) 
and, on the other hand, on the functional dependence At = y4(e(t)). 

With this approach, we obtain generic results which are not asymptotic in time, but valid 
at any time t and for a large class of functional dependence on the uncertainties. 

Though the approach is different, our conclusions are consistent with the cases presented 
in Sect. [2l We extend the observation of Lande that, when adding environmental noise to 
growth rate on the log scale, environmental stochasticity increases the mean multiplicative 
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growth rate to matrix models. As to Tuljapurkar's asympotic approximation, we arrive at 
a different conclusion because his assumptions correspond to a matrix A depending linearly 
on the environmental factors as in ([2]), and our result does not cover this case. 

Our general conclusion is, therefore, that the analytical dependence on environmental 
factors cannot be overlooked when trying to tackle the influence of variability. However, 
as shown in this paper, specific analytical dependence coupled with appropriate notion of 
variability lead to wide generic results, valid for all times and not only asymptotically, and 
requiring no assumptions of stationarity, of normality, of independency, etc. 
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